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Abstract 



Certain perturbation graphs in the calculation of the effects of the medium 
on neutrino scattering in supernova matter have a nonintegrable singularity in 
a physical region. A number of papers have addressed the apparent pathology 
through an ansatz that invokes higher order (rescattering) effects. Taking the 
Gamow- Teller terms as an example, we display an expression for the spin-spin 
CT} . correlation function that determines the cross-sections. It is clear from the form 

that there are no pathologies in the order by order perturbation expansion. 
qq ' Explicit formulae are given for a simple case, leading to an answer that is very 

Q\ ■ different from one given by other authors. 

Q_i' There is a considerable literature 1-6 directed at the time dependent spin fluctua- 

tions in the hot nuclear matter of the postshock supernova environment, with appli- 
cation to the neutrino physics in the medium. The nucleon spin correlation functions 
are of particular importance in these problems because the Gamow- Teller part of the 
weak interaction supplies roughly 3/4 of the rate for neutrino- nucleon processes. The 
spin (and spin-isospin) correlation functions determine the GT reaction rates and are 
known to be strongly affected by the N-N interactions in the medium 7-11 . However, 
there has been disagreement both about the relative importance of the different cor- 
rections, and about the detailed mechanics of calculation. Much of the work referred 
to above attempts to go beyond perturbation theory; but one of the disagreements 
about method is rooted in the properties of perturbation theory. This has to do with 
the nonintegrable perturbative singularities in the time Fourier transform variable, u, 
coming from simple graphs in an approach based on T matrices. The authors of refs 
0-|J advocate the regulation of these singularities by rescattering affects that they 
estimate through a physical rather than an analytical argument. In contrast, I have 
argued that these singularities are regulated, in the calculation of neutrino scattering 
effects, by contributions from "medium dependent wave function renormalization" 12 . 

In the present note we consider the Fourier transform of the time-dependent spin 
correlation function in an interacting system, and present a formal expression that 
is well adapted for the demonstration that the singularities must disappear, order 
by order, in physical results. The end results are the same as those provided by 
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a T matrix approach followed by regularization via medium-dependent wave func- 
tion renormalization, but the latter procedure is no longer required. Our procedure 
is also simpler, both conceptually and in the application at hand, than is the con- 
struction in the imaginary time formalism of the requisite thermal Green's functions, 
followed by the continuations needed to get the real time correlation function. When 
the spin densities are replaced by the complete set of weak current operators, the 
approach laid out below also is well adapted to perturbative consideration of the cor- 
relation functions for all of the currents. It should be especially useful in resolving 
the disagreements that persist in the calculation of medium dependent wave function 
renormalization that has been needed in calculating the effects of the photon bath on 
the weak interaction rates in the early universe, just before nucleosynthesis 13-16 . 

The Fourier transform of the correlation functions for the spin densities of an N 
particle system is given by, 

S ifj {q,w) = Z- 1 J d 3 xdt e _iq - x e <wt 7V[e- /JH <7i(x, t)a 3 (0, 0)}, (1) 

where the <7j(x, t) is the spin density operator in the Heisenberg picture, and Z is the 
partition function. 

To determine the rates of neutrino reactions in the medium this correlation func- 
tion is integrated with respect to a weight function, gi t j(q,u>), 

W = J2j d 3 qduj Si j(ci, w)gy(q, u), (2) 

hi 

where (q, u) is the momentum-energy transfer to the leptons and where the weight 
function is determined kinematically, depending on the reaction under consideration. 
We need to distinguish two cases: 

1. The support of the function, g(q,u) includes a region around u = 0; 

2. The function g(q, to) vanishes for u < or for uj > 0. 

Neutral current v scattering is an example of case #1 and v pair production or 
absorption are examples of case #2. 

To extract terms in the perturbation development of ([I]) we make a transformation 
to an interaction picture based at time t = 0, 

(n(x,t)= U(0,t)a[ I] (^t)U(t,0), (3) 
where U(ti,t2) is the interaction picture time translation operator, 

U(t 2 ,t 1 )=T[exp(-t f 2 dtHj(t))] = 
Jti 

I -i f 2 dt'H^t') - f 2 dt'H^t') f dt'H^t") + (4) 
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Hi(t) is the interaction Hamiltonian in the interaction picture. In (f|) we write the 
familiar perturbation expansion. It defines a continuation for complex (ti,t 2 ), in the 
case in which Hj(t) is an analytic function of t, so that the integrals are independent 
of path. With the continuation so defined, we also have the result, 



e-* H = e 



- pHo U{-ip,0), (5) 



which can be easily verified by observing that the continuation of (Q) provides the 
"r-ordered" operator from finite- temperature equilibrium quantum mechanics 17 , 

U(-i0, 0) = T[exp(- f dTHA-ir))}. (6) 

Jo 

The interaction picture spin-density operator, a\ (x, t) (which is the axial vector cur- 
rent for nonrelativistic nucleons), can be written in terms of the Pauli matrices for 
the individual particles, af^ 

N 

oj(x, = ^ 5(x - x a - mT^at + c a )af\ (7) 

a=l 

where x Q is the Schrodinger position operator and p a the momentum operator for 
the a'th particle and the c Q 's are irrelevant constants. We can use and (|5|) to 
write (1) as 

S id (q,u) = Z- 1 J dxdt e- iq V w 'Tr[e-^°[/(-i/5,t)^ /) (x, t) U{t, 0)a { /\o, 0) 

where we have used the composition properties of the U operator, U(—i/3, 0)U(0, t) = 
('( '>'•/)• 

Eq. (^) is our basic result. The interaction Hamiltonian enters only through the 
U operators, and it remains to expand them in the power series (ffl) and take the 
matrix elements. We note some features before giving an example. 

A) Having based the interaction picture at t=0, we will generate more terms in 
a given perturbation theory time integral than one encounters in constructing the 
correlation function from on-shell T matrix elements weighted by statistical factors. 
This is because the T matrix elements are calculated in an interaction picture based 
at t — — oo, with adiabatic cutting off, and in a typical time integration one therefore 
gets a contribution only from the upper limit, whereas in our time integrals in the 
expansion of matrix elements of U there are finite contributions from both limits. 
The cost of the T matrix approach, is considerable, however, in requirinq a "medium 
dependent wave-function renormalization" 12 , that must be done carefully. The need 
for this procedure is a consequence of turning the interaction off asymptotically when 
the interacting particles do not become separated at long times. This is analogous 
to the case of a field theory in a vacuum, where the self interactions, which do not 
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disappear at large times, lead to a wave function renormalizations for all participating 
particles, with the difference that now it is the interactions with the medium and not 
the self- interactions that must be taken into account. 

As remarked above, if one follows a T matrix approach, but without including 
this wave function renormalization , and calculates W, for case #1, one encounters 
non-integrable singularities in the lowest- order perturbation correction to the free- 
particle correlation function. In contrast, the basing of the interaction picture at 
t=0, as in ([8]) avoids both any singular terms in the straightforward order by order 
expansion of W, and the conceptual complications of the T matrix approach as well. 

B) The operator U(—i/3,t) in (|S|) combines part of what is in T matrix language 
the "scattering" correction (the t to part) with the correction to the statistical factor 
(the to —if3 part). We can now do the integrations for the combined factor directly, 
saving in the number of terms generated. Thus it seems best on all counts to use 
the form,(||), with the interaction picture based at t = 0, as the basis for calculating 
correlation functions perturbatively. 

As an illustration, we take the model of Raffelt and Sigl 1 , where only the auto- 
correlation function of a single particle that moves in one dimension is considered 
and where the interaction is through a fixed spin-dependent potential. In ref. [0, 
this potential, designed to simulate the rest of the medium, is one produced by a set 
classical moments of random direction, situated at random points, x a , each with an 
(interaction picture) interaction with the particle, 

< s,p\Hi a) (t)\s',p' >= v(p-p') jl a ■ a SyS , e <P-P')-*« e i{E P -E p ,)t^ (Q) 

where fi is the classical moment and the states of the particle are labelled by a single 
momentum and spin, p, s, and E,E' are the energies of the states labelled by p,p'. 
For our demonstration, it suffices to take a single one of these sources. At the end 
of the calculation the correlation function will get averaged over all positions of this 
source, so that the implicit use of translation invariance in ([]]), in using the points, x 
and 0, is allowed. In ref.|l|] the authors find that their results are quite independent 
of the momentum transfer, q, over the range of q defined by the weight function, 
gij(q,uj). For our tests we set q = in the correlation function; the formulae below 
would be slightly more complex if we kept the q dependence, and the conclusions the 
same. Then, using (|7|) with N=l, we perform the x integral in ([8]) obtaining for the 
second order contribution, 

Sij(q = 0,u) = Z- 1 J dt e^'Tr[e- /3Ho f/(-z/3,t)(T i f/(t,0) ( T j ], (10) 

where Oi is simply the Pauli matrix for the single dynamical particle. To evaluate 
(|l0l) to second order in the strength parameter, fi, we write the expansion (f|) as 
U = + + and perform the time integrals, obtaining, 
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W = Z-\2tt)- 1 J dqdpduj gi j{q, u) J dte iuJt < p, s\e~P H ° [aiU {2 \t, 0)o>- + 

i,j s=l,2 

U {2 \-ip,t)<j i( J 3 + U {1) (-il3, t)aiU {1) (t, OyJ \p,s > . 



[ID 



Using the interaction Hamiltonian 
integrals we obtain 



in the terms of (01) and performing the time 



W 



* ^5{u-Ei+E)+[2P(E x -E)-l§\8(u) 



(Ei-E) 2 



7T 



12 f o!pdpie-^ g l 
" J Kf>i-p)| 2 



f gi (Ei -E)+[2f3(E 1 -E) -f]g 1 (0) 



{E\-E) 2 



(12) 



where <7i(w) = Sii J dqgij(q,uj), and where we have averaged over the directions of 

^. In case #1, for which the prototype is neutrino scattering, we can now see from 
(O) both how it may appear that there is singular behavior that needs regulation, 
and why that is not the case. If in the middle step of flI2| ) we do the p\ integration 
before the u integration, and take just the first term in the square bracket, we do 
indeed get a uj~ 2 singularity, and it is not removed by any other uo~ 2 singularity. This 
first term is the whole of the "perturbation" result quoted in ref. |] and essentially 
that of similar expressions in refs. @-[[|. But if we look at the results of integrating 
over uj first, in the last line of (|P2"D , where the singularity of this first term turns 
into (E — -Ei) -2 , we see how the other terms remove both the second and first order 
singularities, leaving a regular remainder. The physics is simple; in the middle step 
of ( |T2"D the 8{oj — Ei + E) term has to do with the medium being left excited by 
the transfer of energy to, whereas the cancelling term with 5(u) is from the process 
in which the medium is left unexcited. This is the term that is introduced by the 
"medium dependent wave function renormalization" in ref. We see that in the 
present more systematic approach it appears from a completely mechanical calculation 
of the correlation function. Thus we conclude that the replacement by a Lorentzian 
form, uj" 2 — > {uj 2 + A 2 )" 1 , as in refs. [p]] -[|], is neither necessary nor justified. 

Turning to case #2, for which the prototype is neutrino pair production, we 
observe that, to make sense, the weighting function g(u>) must vanish at u = 0; in 
the neutrino pair production problem it vanishes as uj 3 . The 5(u) term contributes 
nothing to this process and the uj~ 2 singularity is harmless. Therefore there is again 
no reason to regulate the u~ 2 terms in the correlation function oj~ 2 . 

To make this point more sharply it may be useful to redescribe the problem that 
is being addressed here. We started with the question of how the interactions of 
nucleons in the medium affect the neutrino scattering rate. Suppose we ask instead 
how the interaction with a neutrino would affect the scattering of a nucleon within 
the matter. Our amplitudes would have a lo~ x singularity, and the rate corrections 
an u~ 2 singularity, the latter again being removed, in an inclusive calculation, by a 
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non-analytic term, S(u). Notwithstanding this order-by-order removal in the inclusive 
rate, which has to do with the interference of connected and disconnected T matrix 
graphs in the construction of the correlation function, let's look at the analytic (and 
totally connected) part of the amplitude by itself. The singularity comes from the 
graphs in which the weak current, which transmits energy u, attaches to an external 
nucleon line in the nuclear process. It is a simple pole singularity in the amplitude, no 
matter how complex an N-N amplitude it attaches to, and the same holds true if we 
consider many-nucleon amplitudes. The point is that in every order in the nucleon- 
nucleon interaction, the singularity is still a simple pole. If the sum of graphs were to 
move the position of the singularity, as in the case of propagator modifications, or to 
exponentiate the singularity, as in the infrared problem, we would have found higher 
order singularites in the higher order graphs. We conclude again that the regulation 
proposed in refs. [Q-Q does not occur. This is not to say that multinucleon effects 
cannot produce interesting modifications of the exact correlation functions for small 
uj. However, we are skeptical that the mock-up of many body effects used in ref.[|IJ, 
in which all but one nucleon is replaced by a set of fixed potentials with randomized 
positions for the sources, can be illuminating in discussing this behavior. 

A numerical nonperturbative calculation is done in ref.[jl[ as a function of u, but 
only for u> ^ 0. The numerical calculation cannot be extended to uj = to produce 
the 8(co) terms that are required to cancel the singularity in case #1, so that the 
calculation itself does not address the issues raised in the present paper. However, 
we can compare with our result ([12]) the result, quoted in |1], for the regulated 
perturbation result in the model. We do this for case #1, taking neutrino scattering 
kinematics, and in one dimension. The ansatz of refs [|lj-|| leads to the replacement 
of (0) by 

w es = z-v-y f dpipie ~Z\\, f3l{E \; E \l d3) 



v( Pl -p) | 2 1 \-(E 1 -E) 2 + I^- 
Using parameters from ref.lJ we take v(p) = vq exp[— (bp) 2 /A], where 6=1F, and 
T = .5(T/lMeV) 1//2 MeV. For gi(u) we take the function that comes from the integrals 
over the initial and final neutrino phase spaces, weighted with the appropriate Fermi 
factors and energy conserving delta function, of the leptonic part of the squared 
Gamow- Teller matrix element. We do this in a one-dimensional world, for consistency. 
Since we are interested only in comparing flTj) and flT3|), we can normalize to <7i(0) = 1. 
We obtain, 

9l (u) = 29(u)(l - e-^)- 1 [log(2) - log(l + e"^)] + 

+20(-w)(l - e^)- 1 [-log(2) + log(l + e^)}. (14) 

We have calculated the ratio of the rate corrections given in (|T2J) and (|T^) for 
a range of temperatures from 2.5 MeV-100 MeV. W is negative and W reg positive 
for all temperatures, with the ratio W/W re9 ranging from —.70 at T=2.5 MeV to 
— .15 at T=100 MeV. Note that these correction terms are proportional to the same 
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unspecified strength parameter. The perturbation theory is justified, at best, only 
when this strength is small enough so that the perturbation is smaller than the 
unperturbed rate; we are not predicting negative rates. 

A similar reduction in rate coming from interactions with the surroundings has 
been found in a related three-dimensional model 12 . Reverting to a T matrix descrip- 
tion, the second order contribution has two parts. The first reflects a change in state 
of the medium and is positive, since it is, loosely speaking, connected to the opening 
of another channel. The second reflects a second order correction to the wave-function 
of the nucleons in the medium, and leaves the medium unexcited, so that there can 
be an interference between zeroth and second order amplitudes. This term is negative 
and outweighs the first term. 

These comments do not apply to the neutrino pair production rate; in this case 
there is no zeroth order term, the leading second order term is positive, and the 
discrepancy between the results of our method and the results of refs.||-[||]is simply 
the change introduced by the presence, in the latter, of the arbitrary regulator, A. 

Finally, we note that the model defined by the single particle interacting with 
classical magnetic moments, with form factors, does not simulate important, and 
calculated, effects of true many-body interactions on current correlation functions. 
The ring graphs, which involve at least two dynamical particles, are not included in 
any sense in the calculations of ref . CJ . The need for the ring sum is not an arbitrary 



singling out of a special set of graphs; as laid out in ref. [11], among other things, the 



ring sum embodies the long-wavelength limits of the four correlation functions that 
determine the neutrino scattering process. The ring corrections lead to enhancement 
or suppression, depending on the sign of the force, if we are doing actual perturbation 
theory with a potential; or depending on the signs of Fermi liquid parameters if we 
follow the approach of [[DJ . This is in contrast to the perturbative terms produced in 
the above model, which do not depend on the sign of the interaction. 

Thanks are due to Adam Burrows for a number of discussions and to Lowell Brown 
for instruction in thermal quantum mechanics. 



7 



References 

[1] G. Raffelt, G. Sigl, |hep-ph 9808476 

[2] G. Raffelt, T. Strobel, Phys.Rev. D55, 523 (1997) 

[3] G. Raffelt, D. Seckel, G. Sigl, Phys.Rev. D54, 2784 (1996) 

[4] H. T. Janka, W. Keil, G. Raffelt, D. Seckel, Phys.Rev. Lett., 2621 (1996) 

[5] G. Raffelt and D. Seckel, Phys. Rev. D52, 1780 (1995) 

[6] G. Raffelt and D. Seckel, Phys. Rev Lett. 67, 2605 (1991) 

[7] N. Iwamoto and C. J. Pethick, Phys. Rev. D25, 313 (1982) 

[8] R.F. Sawyer, Phys. Rev. C40, 865 (1989) 

[9] S. Reddy, M. Prakash, and J.M.Lattimer, Phys.Rev D58,l (1997) 
[10] S. Reddy, M. Prakash, J.M.Lattimer, and J.A.Pons, |astro-ph/98TT294 
[11] A. Burrows and R.F. Sawyer, Phys Rev C58, 554 (1998) 
[12] R. F. Sawyer, Phys. Rev. Lett. 75, 2260 1995 
[13] D. A. Dicus et al, Phys. Rev. D26, 2694 (1982) 
[14] J. L. Cambier et al, Nucl. Phys. B209, 372 (1982) 
[15] R. F. Sawyer, Phys. Rev. D53, 4232 (1996) 
[16] S. Esposito et al, |astro-ph/980"81~96: 

[17] A. L. Fetter and J. D.Walecka , Quantum Theory of Many Particle Systems (New 
York, McGraw Hill) 1971, eqn 24.9 

[18] A. Burrows and R.F. Sawyer |astro-ph / 9804264 



8 



